k 3 superfluid spectrum of highly curved interacting quantum vortices 
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Presented is a prediction, based on the Frenet-Serret differential geometry of space curves, that the 
wave number dependence of the average kinetic energy per unit length of two mutually interacting 
highly curved quantum vortex scales as k~ 3 . The interacting quantum vortices are helical in shape, 
supporting circularly polarized counter-propagating waves, with arbitrary curvature and torsion. 
This power-law spectrum agrees with the high-fc spectrum found in precise quantum simulations of 
turbulent superfluidity with tangle of highly curved and excited quantum vortices. 
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Superfluid turbulence is an intriguing low-temperature 
phenomenon with power-law energy cascades that are un- 
dergoing active investigation. In this Letter we consider 
the origin of a fc~ 3 power-law in the kinetic energy spec- 
trum at high-fc wave numbers (<J the inverse coherence 
length) associated with highly curved quantum vortices. 
The simplest theory for a superfluid condensate in the 
zero-temperature limit is 



ihip*d t if+ - — (VV 
2m 



) • V<p + fi <p*ip - ip*V H <p, 



(1) 

where <p(x) is a complex scalar field for the Bose-Einstein 
condensate (BEC), V H is a local self-consistent Hartree 
potential, and p, is the chemical potential. Minimizing 
the action J d A x £ BEC leads to the Euler-Lagrange equa- 
tion, a nonlinear Schroedinger equation when Vh(M 2 ) = 
ig|(/3| 2 , where g is the real- valued coupling strength of 
the nonlinear interaction, known as the Gross-Pitaevskii 



equation (GPE) 1 ' 2 ihdttp 



It captures complex vortex-vortex interplay (nucleation, 
emission and absorption of vortex rings, and reconnec- 
tion) as well as quantum Kelvin wave mode (kelvon) pro- 
duction on the vortices and sound modes (phonons) that 
escape into the bulk region of the quantum fluid. 

To determine a steady-state solution of the quantum 
vortex in a superfluid, with a BEC wave function de- 
noted by </? v , one solves the time- independent GPE: 

-£ 2 VVv + \ip v \ 2 -l\<p w = with healing length 

£ = h/y/2mfl. A solution for the background condensate 
wave function of a single rectilinear quantum vortex (with 
vorticity along z) is found by separation of variables in 
polar coordinates. Inserting ip^(r, z) — tfi v (r)Z v (z)e 
into time- independent GPE with g/p = £ 2 gives the fol- 
lowing equations with a separation constant k 2 : 



d 2 Z v (z) 
dx 2 

d 2 </> v (r) 1 d(f> v (r) n 



kjZ v (z) = 0, 



(2a) 



dr 2 



dr 



(2b) 

fcjj. Equation ([2a|) admits sinusoidal so- 
lutions and (|2bp can be solved for integer winding num- 



ber n. For the simplest n = 1 case, the Pade approx- 

imant v (r) = y/a 3 $4+Zhi28+ i solves the timc " 
independent GPE with errors at O [(r-y/a) 7 ] . Notice that 
<j) v (r) — » y/a and r — > oo, and thus the nonlinear term in 
(|2bj) vanishes in the bulk. 

The rectilinear quantum vortex solution of the GPE 
is ip,{x) = 4> v {r)[Z+e l ^\\ t+k \\ z ^ + Z-e^w 1 -^} , with 

h 2 k 2 

a parabolic dispersion relation fbu\\ — -o^r- For a vor- 
tex line along the £-direction in cylindrical coordinates 
with unit winding number, the irrotational part of the 
superfluid velocity has a divergent (perpendicular) part 
characteristic of inviscid flow and an advective (parallel) 
part characteristic of rigid translation 

h , s h ' hkn 

v = —V{d±k»z) = 1?± — ±z. (3) 

m mr m 

From this velocity field, wc know that the circulation is 
quantized, k = § dl ■ v = From Stokes' theorem, we 
have f dS ■ V x v = f dS ■ u> = — . For a rectilinear z- 

J J m 

directed quantum vortex, the real part of the vorticity is 
pinned at the vortex center u) — ^5^ 2S> (r)z. For a vortex 
filament of any shape, say a curve C, the velocity field in 
general may be written as 



v(r) 



h f ds' x (r — s' 
2m J c \i 



a f\3 



(4) 



where ds' is the differential length along the vortex fila- 
ment, s' is the parametrization of C, r is the field point, 
and ds' is the differential line element at the vortex center 
and parallel to the vorticity^. The Biot-Savart formula 
(|4]) reduces to v = r dh/(mr) for the case of an infinite 
rectilinear quantum vortex positioned along the center of 
a cylindrical coordinate system with s' = z, ds' = dzz, 
and ds' x (r — s') = -ddz \z\/\r — s'\. 

Using the Madelung transformation 4 <p = y/pe lS ^ h , 
([T]) can be written in terms of the conjugate fluid vari- 
ables (the action S and number density p) as follows: 

+ It(^) 2 -m + ^(p) . 



-pd t S - p 



(vg) 2 

2m 



The quantity in square brackets is identified with a 
Hamiltonian energy functional. Bohm originally made 



2 



this identification while considering quantum flow in 
a spatially-dependent linear external potential V(x)^. 
For a BEC, the nonlinear internal potential energy is 
E int (p) = V H (p) — ijl. A semiclassical energy functional 



is thus H B 



(VSf 



h 

2m 



(if) + E Up)- The average 



energy is a statistical volume integral E — J d 3 x p H BEC 
with p taken as the effective probability distribution 
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E is a con- 



served quantity. Here we are interested in the first term 

^kin = / ^ x P ^2m = / d 3 x p(x) (\mv(x) 2 ) , where the 
de Broglie relation mv = VS is used. The average clas- 
sical kinetic energy per unit length L of a single linear 

quantum vortex is £ = E^ in / L: 



mp 



2 47rm 



— = mp —\og—, 
r 47r r c 

(5) 

where k = h/m is the quantum of circulation and p 
is the constant background number density of the con- 
densate, r\> is a regularizing parameter associated with 
the size of the vessel containing a single quantum vortex 
and r c is an effective cutoff parameter to the divergent 
angular velocity field. 1 

Consider the original treatment by Fetter for nearly- 
rectilinear vortices^. The initial (single vortex) equi- 
librium states are parallel linear filaments with an arc 
length parametrization given by the following vectors 
Ri = (n,Zi) and R 2 °^ = (r 2 ,z 2 ). The consequent de- 
formed state due to the mutual interaction of the vortices 
is parametrized by 

R\ = (n + u 1 (z 1 ),z 1 ) R 2 = (r 2 + u 2 (z 2 ),z 2 ), 

(6) 

where u± and u 2 are treated as small amplitude pertur- 
bations in the radial directions with respect to the initial 
unperturbed filamentary lines. Each vortex filament is 
a stretched helix, approximating a nearly straight line 
parallel to the £-axis. 

The first step is to calculate the fluctuation in position 
of a vortex element at the first vortex due to the presence 
of the second vortex. The fluid velocity at the first vortex 
located at R\ caused by the second vortex at R 2 is given 
by the Biot-Savart law (|4]) 



ds 2 x (R 1 - R 2 ) 
\Ri-R2\ 3 ■ 



(7) 



r c may be chosen so that the integral for £, with a cutoff that 
avoids the singularity at the origin, is equivalent to the origi- 
nal nonsingular integral with no cutoff. So, r c is technically not 
a cutoff parameter per se because the original integral is non- 
singular. Instead, r c is merely a matching parameter useful for 
replacing a nonsingular but difficult integrand with an analyt- 
ically simpler one. All the expressions of the average classical 
kinetic energy per unit length in {3} are equivalent and apply to 
a nearly straight (high-tension and low-curvature) vortex. 



With i = 1,2 denoting the vortices, the differential arc 



length is dsi 



dzi (z + du ^ Zi ^ • Then, accord- 



ing to ([7]), the fluctuation of the position of the vortex 



element originally at R^' is 

du 



dui inj k f 1 - 

= — / dz 2 - 

dt 4?r / 



S ) x 



£zi2 + u\ 2 ) 



1-12 + ZZi 2 + Ui2 



(8a) 



where r i2 = 1*1 — r 2 , Zi 2 = Z\ — z 2 , and u± 2 = Ui — u 2 . 
Also, defining R^ = T\ 2 + ZZ\ 2 and making use of the 

Taylor expansion = ^ 

we can write a leading order expansion of the position 
fluctuation (velocity of the first vortex element at z\ due 
to the presence of the second vortex) 



3u 12 R 1 



dui 
~dT 



a] k _ 
An' 



dzn 



rr 12 + u 12 - z 12 



du-2 

dz 2 



L 12 

3r i2 ■ u 1 2 

" \R?2 ] \ 5 



-ri2 



(8b) 



Making the analogy to the mutual inductance of two 
line currents, the Neumann formula can be used to cal- 
culate the interaction energy of two vortices 



E 12 = mp 



,o v(x) 2 0} mpnl 
2 4tt 



C-2 



dsi ■ ds 2 

\Ri-R2\ 



(9) 

The next step is to develop an expansion for the interac- 
2 (i+±n. .±12.+— \ 



tion energy E\ 2 



\B^+u 12 \ 



in I Jdz!dz 2 - 
where the cross-terms vanish because, in the reference 
frame at the original center of the unperturbed ith vortex 
line along z, the motion of the perturbed filament is in 
the polar direction f^/I^H ~ ±#- Employing Taylor's 
theorem, the denominator is expanded to second order 

1 _ 1 T-12-M12 J_ 3(ri 2 -Ml2) 2 _ (ui 2 ) 2 

111 , 



\R[° 2 ) +u 12 \ 
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so in turn the interaction energy expansion becomes 

1 r 12 ■ u 12 - \{u 12 f 



E\ 2 = 



mpK 
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z\ 2 



dz\dz 2 
du 2 



|< 2 )|3 dz 2 



Ml 



\{r 12 ■ u 12 ) 



(10) 



where in (fT0|) the term obtained by integrated by parts 

1 z -"-12 



was rewritten as 



Hence, it 



^|< 2 »| " |B<° 2 >|3 - |< 2 , 

is straightforward to calculate the variation of the mutual 
interaction energy with respect to a fluctuation at z\ of 
the center of the first vortex line 



SE, 



Sui(zi) 



mpK 2 

47T 



dzy 



r\2 + «12 - ^12 



du 2 

dz 2 



, 3r i2 • «12 

H m — ri2 



(11) 
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Comparing this result with the previous result (18b[) yields 
the useful relation 



mpn - 



dui{z\) 



dt 



Sui(zi) 



(12) 



The mutual interaction energy part of the condensate en- 
ergy arising from a perturbed quantum vortex of length 
L = § c \J dz 2 + du 2 is primarily due to its bending, as- 
suming a sufficient separation distance exists between the 
vortices so that vortex- vortex straining has no low-order 
effect. Therefore, we have 



8E X 



(9) 
(5i 



SL 



d 2 Ui(zi) 



Sui(zi) 



Sui(zi) 



dz{ 



Inserting this result into (|T2|) . yields 



mpKo 



dui(zi) 



dt 



£z x 



d 2 u\(z\) 



dz 2 



(13a) 



(14) 



Having completed our review of Fetter's treatment, let 
us consider a high-curvature generalization of (|14j) . 

The Frenet-Serret formulas of multivariable calculus 
concerning the geometry of curves describe the kinematic 
properties of a particle at position R moving along a con- 
tinuous and differentiable curve C (the particle's trajec- 
tory or world line) embedded in three-dimensional Eu- 
clidean space IR 3 



dt „ dn 

— = nn —— — —nt + to 

ds ds 



where s is the arc length parameter along C, i is the unit 
tangent to C at the point R, h is the unit normal perpen- 
dicular to t at R, b is the unit bi-normal perpendicular 
to both t and n, and n is the curvature and r is the tor- 
sion of C at R. Given a fixed curve C, one constructs a 
local Frenet-Serret frame as follows 



_ R'(s) 



\R'(s) 



i'js) 
\i'(s)\ 



i)Etxn, 



(16) 



where the prime indicates differentiation with respect to 
s. So n points along the direction of the derivative of 
t with respect to the arc length parameter of the curve 
and equating (jT5j) with (fH)|) . the curvature is k = \t'(s)\. 
Therefore, the unit vectors t, n, and b serve as an orthog- 
onal coordinate system centered at R, a local reference 
frame that moves with the particle. For example, spec- 
ifying the points in K 3 with the cylindrical coordinates 
(r, i9, z), if C is a helix with its axis along r = 



R(s) = (rcosfcz,rsinfcz,0)+(0,0, t)kz) (17a) 



u(kz 




(17b) 



with kz = s/ y/x 2 + f} 2 , then the Frenet-Serret frame is 



ri9 + f)£ 



b=-W + t Uis) 



^/t 2 + f) 2 y/v 2 + \) 2 

and the local curvature and torsion of C are 
r _ f) 



I) 2 



I) 2 



(19) 



In the limit [)<t, the curve is a compressed helix where 
each cycle approximates a circle with curvature 1/r. In 
the opposite limit t) ^ r, the curve is a stretched helix 
(approximating a straight line) with torsion 1/f). 

We may reconsider the case of a perturbation to a rec- 
tilinear quantum vortex: a quantum vortex supporting 
a small amplitude circularly-polarized plane Kelvin wave 
counterrotating in a sense opposite to the vorticity direc- 
tion of the unperturbed vortex line. The solution we seek 
is based on the general Frenet-Serret formulas. 2 Since the 
second derivative of the radial perturbation of the quan- 

turn vortex center is u'{ = —k 2 Ui, (fl"4"|) takes the form 
of an undamped Bloch equation 



dui(zi) 
dt 



Ek 2 
mpK 



■ ux(zi) x s 2 . 



(20) 



The radial displacement iti behaves like the magnetiza- 
tion vector of a nuclear spin precessing about a back- 
ground magnetic field along s 2 ~ z. Thus, two segments 
of the mutually interacting quantum vortices behave like 
coupled nuclear spins. 

The helix rotates in time as a sinusoidal perturbation 
Mi ~ ti e l ( fc2 -"') t so (IT41 may be written as 



iuiui(zi) 



£(k)k 2 
mpKo 



Z X Ui(zi). 



(21) 



Then with U\ — (x, y), we have iuj{x, y) 



£(k)k 



or in matrix form 



ek 2 



iuj 



(;) - (s 



The well 



2 To make additional contact with the previous literature, con- 
sider the limit of small curvature with the perturbed vortex is a 
stretched helix. The local frame, which we may choose to fix at 

zi centered on the first quantum vortex, is i = ' 
Vl ' (zi) =-# + ... , and b 



t x ft = — 19 + ■ ■ ■, where the 
The mutual interaction of 



- |t»i'(*i)| 
arc length is parametrized by s z\ 
the vortices causes them to bend into the filamentary shape of 
a rotating helix — a circularly polarized Kelvin wave. Equation 
H14B may be rewritten as 

. tie! s u' „ 



mpKo \u'\ 



In a parametrization with |it'| = 1, this is the local induction 
approximation (LIA) used for quantum turbulence simulations 3 . 
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known solution by setting the determinant to zero, for 
rjj = 1/k, is the Kelvin wave dispersion relation 

g(fc)fc 2 JSJ n Q P 1 

uj k = = — — log—, (22) 

mpK 47r r c k 

valid in the limit where t/f) is a small quantity with C 
nearly a line. Now, with a locally helical-shaped quantum 
vortex parametrized by (|17[) . we need not restrict t/f) to 
be a small parameter. Instead the only generic constraint 
imposed upon the shape of the perturbed quantum vor- 
tex is that it be locally parametrized by (|17al) . We will 
not loose the generality of two interacting quantum vor- 
tices of arbitrary shape by restricting their mutually in- 
teracting segments to helices. 

Since \R'(s)\ = 1, i = R'{s) and in turn jUJ is 



d 2 R(s) (J! t 



-n. 



(23) 



ds 2 t 2 + f) 2 

So, allowing for high curvature, (|16p can be written as 

t 2 + f) 2 . d 2 R(s 



-t x 
t 

= -t x it(s) 
t 



(24a) 
(24b) 



-0 + xz ® ^ + JzJxE(s), (24c) 

We are now considering a pair of parallel vortex segments, 
arbitrary helices equal in magnitude and spin, circling 
around each other whereby each segment moves by the 
influence of the velocity field of the other. Since the veloc- 
ity field of a vortex segment (with a circularly-polarized 
radial perturbation) is anti-parallel to the bi-normal unit 
vector of a helical curve (i.e., v cx — b), this velocity field 
may be generally expressed in the Frenet-Serret frame as 



v = uj {i)xt — xz) 



(25) 



Next, let us determine the radius t and the stretching 
f) of the helix directly in terms of the available physi- 
cal parameters describing the mutually interacting quan- 
tum vortices, such as the separation distance between 
the vortices 2R^ ■ Let us denote the orbital radius at 
R = R[ 2 \ Thus, the magnitude of the polar velocity 



component is v& 



i(2R a 



If v& = lj x R , where 



u> = uj z is the effective orbital "cyclotron" spin vector, 
then the orbital angular frequency of the vortex pair is 

/i — £jL — h 
^° ~ Ro ~ ImKi 



4 *j^ a . The quantum vortex solution 
([3]) with a plane- wave phonon mode counter-propagating 
to u) along the — z direction with wave number k = kii 
is concomitant to the quantum vortex kelvon mode 



|3j H - hk „ 

= 17 Z = U! 

2mR m 



R tf - IRikz 



(26) 



Equating ([26| to (I251) . the helical parameters f) and t 
are thus analytically determined to be f) = R Q and 
t = 2R 2 k. In turn, the curvature and torsion of the helix 

are K = T+4fc% and r = r; i+WKi ■ 0ne observes 
that the maximum curvature K max = k occurs when 
k = 1/(2R ) and the maximum torsion r max = 1/R a 
occurs when k = 0. Both the curvature and torsion van- 
ish in the limit of infinitely separated straight-line vor- 
tices. Inserting (|2"5j) into (|24cl) . one finds a precise self- 
consistent equation governing the dynamics of a quantum 
vortex with large-amplitude helical wave with arbitrary 
wave number k in mutual interaction with another quan- 
tum vortex 



R= -LUai'd 



2R k 



x R(s). 



(27a) 



For analytical continuation to Fetter's treatment in the 
high-tension limit f) 3> t, the fluctuation of (I27a[) obeys 



u 



2R Q k 



(27b) 



Equating (|27b[) to the Bloch equation (f2TJf provides a way 
to determine the wave number dependence £ = £(k): 



£(k) 



mpn uj 
2R k 3 



mpK 
%-kRI 



(28) 



A helical wave triggered by a fc-mode axial phonon 
counter-propagates along the segment undergoing mu- 
tual interaction. Only the inital axial phonon dispersion 
relation lu = f^fc 2 is consonant with the dispersion rela- 
tion for lu k given by (|22|) for a semiclassical Kelvin wave. 



Equation (|28|) is an analytical prediction, based on the 
Frenet-Serret differential geometry of space curves, that 
the helical wave spectrum scales as k~ 3 for highly curved 
quantum vortices. This power-law was found in precise 
quantum simulations of supcrfluid turbulence support- 
ing highly curved vortices^. A k~ 3 spectrum also arises 
from a single rectilinear quantum vorte>£, and so it has 
been recently suggested that this straight vortex power- 
law underlies the high-Zc part of a turbulent superfluid 
spectrum^. The analysis presented here suggests that 
(j28f . not ([5]) originally discovered by Lord Kelvin in 1880 
nor the Fourier spectrum of a single rectilinear quantum 
vortex, may be responsible for the high-/c spectrum be- 
cause of the effect of mutually interacting highly-curved 
quantum vortices characteristic of supcrfluid turbulence. 
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